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Motivated by its potential application to earthquake statistics, we study the exactly self-similar 
branching process introduced recently by Vere-Jones, which extends the ETAS class of conditional 
branching point-processes of triggered seismicity. One of the main ingredient of Vere-Jones' model 
is that the power law distribution of magnitudes m' of daughters of first-generation of a mother of 
magnitude m has two branches m' < m with exponent /3 — d and m! > m with exponent (3 + d, 
where (5 and d are two positive parameters. We investigate the condition and nature of the sub- 
critical, critical and super-critical regime in this and an extended version interpolating smoothly 
between several models. We predict that the distribution of magnitudes of events triggered by a 
mother of magnitude m over all generations has also two branches m' < m with exponent (3 — h and 
m' > m with exponent /3 + h, with h — dy/1 — s, where s is the fraction of triggered events. This 
corresponds to a renormalization of the exponent d into h by the hierarchy of successive generations 
of triggered events. The empirical absence of such two-branched distributions implies, if this model 
is seriously considered, that the earth is close to criticality s ~ 1 so that /3 — h ~ /3 + h ~ @. We 
also find that, for a significant part of the parameter space, the distribution of magnitudes over a 
full catalog summed over an average steady flow of spontaneous sources (immigrants) reproduces 
the distribution of the spontaneous sources and is blind to the exponents /3, d of the distribution of 
triggered events. In summary, the exactly self-similar Vere-Jones model provides an attractive new 
approach to model triggered seismicity, which alleviates delicate questions on the role of magnitude 
cut-offs in other non-self-similar models and provides interesting predictions which could be tested 
with stochastic reconstruction methods of earthquake sequences. 

PACS numbers: 64.60.Ak; 02.50.Ey; 91.30.Dk 



INTRODUCTION 



Stochastic branching processes describe well a multi- 
tude of phenomena 0, from chain reactions in nuclear 
and particle physics, material rupture, fragmentation 
and earthquake processes, to population and biological 
dynamics, epidemics, economic and social cascades and 
so on. Branching processes are also of particular in- 
terest because deep connections have been established 
with critical phenomena. In branching processes, var- 
ious quantities exhibit power law distributions at crit- 
icality. This includes the distributions of cluster sizes, 
of the number of generations before extinction and of 
durations. 

Branching processes become critical when the average 
death rate is exactly compensated by the average growth 
rate. At criticality, branching processes become asymp- 
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totically self-similar, which translates for instance into 
an asymptotic power law tail for the distribution of clus- 
ter sizes. Such scale- invariance is a general characteristic 
of systems at the critical point of a phase transition or of 
a bifurcation. But self-similarity in general holds only 
(i) at criticality and (ii) asymptotically, i.e., at scales 
much larger than the microscopic mesh or elementary 
branch scale. 

Vere-Jones has recently introduced a class of self- 
similar branching processes which is exactly self-similar 
for a broad range of parameters (of non-zero measure), 
that is, far from criticality and for all scales Vere- 
Jones self-similar branching process is derived from a 
class of exactly self-similar random measures which gen- 
eralize the class of stable purely atomic completely ran- 
dom measures. The underlying idea is that a change of 
scale is balanced by a change in mass, making the mea- 
sure, and the branching process, self-similar at all scales 
and off-criticality. This is possible only for branching 
processes with continuous masses or "marks." Then, 
the model has a natural application to describe earth- 
quake triggering, which is the example we adopt in the 
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following to formulate the problem and present our re- 
sults, without loss of generality. 

The concrete example proposed by Vere- Jones is 
the "self-similar ETAS" (epidemic-type aftershock se- 
quence) model, which is a self-similar extension of the 
initial non-self-similar standard ETAS model [(J |20| . 
The main statistical properties of the standard ETAS 
model are reviewed in [j, Q • The ETAS model describes 
the rate A(r, t, m) of events (earthquakes for instance) 
at position r at time t and of mass or mark (magnitude) 
m resulting from spontaneous sources ("immigrants") 
and from all past events. The model is defined in terms 
of the conditional Poisson intensity A(f, t, m) which is a 
function of all past events. The standard ETAS model 
describes earthquake triggering by past ones within the 
framework of branching theory and takes into account 
the interplay between the exponential productivity law 
and the Gutenberg-Richter (GR) law of event sizes. A 
series of papers has shown that the standard ETAS 
model gives a reasonable description of the statistics of 
seismic clustering 6], of foreshocks 171 |8i and aftershocks 
[J HI triggered by other earthquakes, of the empirical 
Bath's law for the largest aftershock of a given sequence 
0, and of the statistics of seismic rates 0. The 
standard ETAS model thus constitutes a powerful null 
hypothesis to test against other models 0. 

The standard ETAS model is perhaps the simplest of 
a much larger class of models embodying the physics 
of triggered seismicity. Its simplicity results from 
its assumption of complete decoupling between the 
Gutenberg-Richter distribution of event sizes, produc- 
tivity law, time and space interaction kernels. Al- 
ready in 1988, Ogata proposed an extension allowing 
for magnitude-dependent time and space kernels and 
taking into account spatial anisotropic effects 6]. Us- 
ing statistical likelihood methods, this extension has re- 
cently been shown to be superior to the standard ETAS 
model to account for empirical earthquake clustering 
[l2|. Ogata et al. have also extended ETAS to charac- 
terize regional features of seismic activity in and around 
Japan, by allowing the parameter values to vary from 
place to place 01 • ^ n the same vein, we have shown that 
the empirical distribution of seismic rate in California 
can be described adequately by ETAS model only when 
taking into account a strongly non-uniform fractal distri- 
bution of spontaneous sources . Ouillon and Sornette 
have considered a multifractal model of triggered seis- 
micity which predicts magnitude-dependent exponents 
for the Omori law of aftershock decay rates, in good 
agreement with empirical data 0, ^| . Unfortunately, 
few of these extensions have allowed yet for a full theo- 
retical understanding of the stationary properties of the 
resulting catalogs of triggered events. 

Here, we present a detailed theoretical understanding 
of some statistical properties of the self-similar Vere- 



Jones model, which extends the standard ETAS model 
by allowing for some dependence between the random 
magnitudes of progenies and ancestors. In an effort to 
clearly formulate some universal key features of the self- 
similar Vere- Jones model, we shall consider mostly the 
statistics of total numbers of aftershocks triggered by 
spontaneous sources and their distribution. 

The organization of the paper is as follows. Section 
2 defines Vere- Jones' model together with a general- 
ized version that allows us to interpolate smoothly be- 
tween the standard ETAS model and Vere- Jones' model. 
Section 3 presents the general theoretical treatment in 
terms of generating probability functions (GPF) which 
predicts in particular the exact form of the magni- 
tude distributions. Section 4 explores the phase di- 
agram of the model by identifying the conditions for 
sub-criticality, criticality and super-criticality. Section 
5 derives the distributions of earthquake magnitudes 
for Vere-Jone's self-similar model. Section 6 presents 
a summary of our predictions, a discussion of their con- 
sequences and their possible empirical tests. 



II. VERE-JONES SELF-SIMILAR ETAS 
MODEL 

A. Definition of the "generalized Vere-Jones 
ETAS model" 

Here, we integrate over space and time to focus on 
global properties, such as the total number of events 
triggered by a given event. Let us consider first a model 
more general than the Vere-Jones "self-similar ETAS" 
model, which we call the "generalized Vere-Jones ETAS 
model." In this model, an event of magnitude m may 
trigger another event of magnitude ml according to the 
Poisson rate X(m,m'), so that the total rate of produc- 
tion of events of magnitude ml is 

A(m') = ^A(m,m') (1) 

m 

where 

A(m,m') = Ke «m-Pm'-d\m- m '\ > ^ 

gives the number of events of magnitude ml triggered 
directly (on first generation) by an event of magnitude 
m. The space and time dependence of the mother (or 
triggering event) of magnitude m and daughter (or trig- 
gered event) of magnitude ml have disappeared from the 
expression of A(m, ml) due to the integration over space 
and time, so that we are concerned only here with total 
numbers in a fixed space-time window. The exponential 
in (0) contains three contributions: 
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1. e om describes the exponentially growing produc- 
tivity of a source as a function of its magnitude 
to 0, 0| (in other words, e om is proportional to 
the average number of first-generation progenies 
of the source to). 

2. e - ' 3 ™ is the so-called Gutenberg-Richter distri- 
bution of the magnitudes of the first-generation 
triggered events, usually improperly referred to in 
the seismological literature as "aftershocks" and 
which is now progressively replaced by the con- 
cept of triggered earthquakes. With this caveat 
in mind, we shall sometimes use the term "after- 
shocks" in a broader sense than used in the stan- 
dard literature to refer to triggered events of arbi- 
trary magnitudes without the condition that the 
aftershock magnitudes have to be smaller that the 
source magnitude. 

3. The new term e - d \ m ~ m 1^ compared with previ- 
ous models of the ETAS class, describes a bind- 
ing or localization of the magnitude to' of trig- 
gered events in the neighborhood of the ances- 
tor's magnitude to. In other words, this term 
means that daughters' magnitudes keep a memory 
of the size of their mothers: large (small) daugh- 
ters come more probably from large (small) moth- 
ers for d > 0. This could be associated with the 
fact that triggered events occur on patches of the 
mother's fault rupture with large residual stresses 
(for the e - d ( m ~ m ) branch with m' < to) and on 
faults branching from the mother's rupture (for 
the branch e d ( m - m ) with mf > m). 

The branching model is such that any event can trig- 
ger other events according to the rate A(to, to') given 
by Thus, a given event may give daughters of first- 
generation, which can themselves trigger other events 
and so on, giving rise to a triggering cascade. The model 
also considers the existence of spontaneous sources (im- 
migrants), seeding the branching cascades. 

B. The ETAS model: d = 

For d = 0, i.e., when the size of daughters is indepen- 
dent (has no memory) of the size of the mother, model 
© recovers the standard ETAS model 0,0,113]: 

A(m, m') = n'e a{m ~ ma ^ x (3 e -^ m '- m ^ , (3) 

where the first factor in the r.h.s is the so-called ex- 
ponential productivity law and the second one is the 
Gutenberg-Richter law of first generation events magni- 
tudes which is normalized (J^ dm' (3e~^ m _m °) = 1). 
A minimum event size toq is necessary to make the 



ETAS model well-defined [UHi]. Indeed, it is not pos- 
sible to make the model convergent for mo — > -co (a 
— oo magnitude corresponds to a vanishing energy, since 
the magnitude is proportional to the logarithm of the 
energy) , and an ultra-violet cut-off is necessary pH |22| . 

The constant factor k' is derived from the coefficient 
k of the generalized Vere- Jones ETAS model as 



The ETAS model is critical when its average branching 
ratio (equal to the number of progenies averaged over 
all mothers' magnitudes) 

n= dmn'e a{m - mo) x ^-/JCm-mo) = (5) 

is unity. The case n < 1 (respectively n > 1) corre- 
sponds to the sub-critical (respectively super-critical) 
regime. The condition a < (3 is needed to make 
the model convergent, as both the borderline a = 
and the regime a > (3 leads to finite-time singulari- 
ties with stochastic times |2i|. Alternatively, conver- 
gence and stationary is obtained by adding an upper 
magnitude cut-off m max , associated with the empirical 
bending down of the Gutenberg-Richter law for magni- 
tudes larger than about 8, suggesting that m max w 8 — 9 
lljl HE US ■ We shall not consider the influence of this 
upper cut-off whose impact is rarely felt only at time 
scales and for earthquake numbers so large to ensure 
that the largest earthquakes are sampled. 

The ETAS model, unlike Vere- Jones model, lacks self- 
similarity due to both the non-generic condition n = 1 
for criticality and the existence of a minimum event 
magnitude mo which introduces a characteristic mag- 
nitude scale. 



C. The self-similar Vere-Jones model: d > and 

a = /3 

Vere-Jones self-similar ETAS model corresponds to 
a = j3 and d > in (gj: 

X(m,m') = A(m-m') , A(m) = ^"^M . (6) 

The condition a — [3 expresses a balance between the ex- 
ponential small probability of finding a large mother and 
its exponentially large productivity. When a — (3 (and 
in absence of the effect of d), each event magnitude range 
[to, to + dm] contributes equally to other event trigger- 
ing, and, in particular, small events are as important 
to event triggering as are larger ones, a propert y w hich 
seems to be approximately true for earthquakes ^18, 19]. 
In the standard ETAS model, the case a — (3 is not pos- 
sible without the introduction of both an ultra-violet 
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cut-off too and an additional infra-red cut-off TOp^x trun- 
cating the Gutenberg- Richter distribution pH l25l l2a | , 
which are necessary in order to obtain non-diverging se- 
quences. In Vere- Jones self-similar ETAS model, the 
condition a = is made possible without cut-offs by 
the introduction of the parameter d > 0. Physically, 
expression © (together with the standard Gutenberg- 
Richter law) can be interpreted by the existence of two 
branches for the Gutenberg-Richter distribution of the 
~ e 0™ events triggered by a given mother of magnitude 



to : 



1. daughters with magnitude 



< 



have 



their magnitudes to' distributed according to 



-Q3-d)m while 
2. daughters with magnitude 



> 



have 



their magnitudes to' distributed according to ~ 

e -(l3+d)n/_ 

The former (respectively latter) distribution branch of 
daughter magnitudes ensures that there are less small 
(respectively large) daughters than with d = 0. This is 
the origin, as we will make clear below quantitatively, 
of the possibility to avoid any cut-off too or TO max and 
still obtain a convergent model. This is a first reason 
why Vere- Jones' model © is self-similar. As we shall 
see below, the other reason is that there is a finite range 
of parameters for which Vere-Jones is effectively critical. 



III. GENERAL THEORETICAL 
FORMULATION 

A. From Bernoulli to Poisson Statistics 

One of the main ingredients of both the standard 
ETAS and Vere-Jones' models of aftershocks trigger- 
ing is the Poissonian statistics governing the number of 
first generation aftershocks triggered by some sponta- 
neous source of given magnitude m. In order to obtain 
a deeper insight into the origin of the underlying Pois- 
sonian law, we start with the more general Bernoulli ap- 
proach to the description of aftershocks triggering statis- 
tics. 

Consider the Bernoulli version of a generalized model 
in which each spontaneous source of magnitude to has 
p independent "possibilities" to trigger some aftershock. 
The probability that an aftershock of magnitude in the 
interval [to', to' + dm'] is actually triggered along one of 
these p paths is denoted D(m,m' ,p)dm' . Keeping the 
approach general, we allow D(m,m' ,p) to depend on 
the source magnitude to, the number of channels p and 
the triggered magnitude to' in an arbitrary way, in all 
what follows in this section. The next section will then 
apply our formalism to the specifications J5J and ©, 



using the standard ETAS version © as a reference and 
point of comparison. 

The generating probability function (GPF) of the ran- 
dom number of first generation aftershocks, obeying to 
the above Bernoulli statistics, is given by 



Qi(z;m,p) = 



1+ dm'D(m,m',p)(z-l) 



m„ 



(7) 



Here, too is the smallest magnitude of possible triggered 
aftershocks, which will be pushed to — oo in the self- 
similar Vere-Jones version. 

An essential assumption of the models studied here is 
that triggered events of first generation act themselves 
as sources which trigger their own aftershocks according 
to the same laws. Let us call Q(z; m',p) the GPF of the 
random number of events of all generations triggered by 
a first-generation daughter with magnitude to'. Then, 
the GPF of the number of aftershocks triggered over 
all generations by a given mainshock of magnitude to is 
solution of 



0(z;to,p)= 1+ I dm' D(m,m' \p)(zQ(z;m' ,p) — 1) 

(8)" 

obtained from (J7J by replacing z by zQ(z; m',p) to ex- 
press that each branch has the same statistical cascade 
properties. 

Consider the limiting intensity of aftershocks trigger- 
ing 



lim pD(m,m! ,p) = A(to,to') 

p — >oo 

and its associated Poissonian GPF limit 



0(z;to) = lim Q(z;m,p) . 

p — >oo 



(9) 



(10) 



Expression (JHJ) leads to the nonlinear integral equation 



0(z; to) = exp 



dm' A(to, m')(z8(z; m') — 1) 



(11) 

It will be useful in the sequel to study the statistics of 
those triggered events with magnitude larger than some 
threshold [A, whose corresponding GPF for their num- 
bers is denoted 0(z; to, fi). The equation for 0(z; to, fi) 
is obtained from 111(1 by replacing 0(z; to) by Q(z; to, h) 
in both the l.h.s. and r.h.s. and by replacing z in the 
r.h.s. by 



H{ji - m') + zH(m' - fi) = 



1 if to' < /! 
z if to' > /i 



(12) 



where H(x) is unit step function, equal to 1 if x > 
and otherwise. Replacing z by l|12fl just means that 
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only those first generation aftershocks, whose magni- 
tudes are larger than /Lt, are counted. Their subsequent 
cascade above the magnitude level \i is accounted for by 
replacing Q(z;m') by Q(z; m', //). This leads to 



0(z; m, fi) — exp [ / dm'X(m, m')(Q(z; mf, fjf) — 1) + 

J mo 
/>oo 

(z — 1) / dm'X(m 7 m')Q(z; mf, pj\ . 



B. Distribution of events magnitudes 



One of the simplest and most informative statisti- 
cal characteristics of branching processes is the average 
number (R) (m, fjf) of events of magnitude above /i trig- 
gered by some spontaneous source of magnitude m: 



(R)(m,fi) = 



dO(z; m, jj) 



dz 



(14) 



Using expression (|13|) in l|14|). we find that (R)(m,fi) 
satisfies to linear integral equation 

/>oo 

(R)(rn,fi)= / X(m,m')(R)(m' , fx)dmf + (Ri)(m, fx) , 

J mo 

(15) 

where 



{Ri)(m,n)= / X(m,m')dm' 



(16) 



is the average number of first generation aftershocks 
with magnitudes larger than /i. 

We assume that the spontaneous sources constitute a 
stationary point process with Poisson statistics, with an 
average number of spontaneous sources during a time in- 
terval t equal to lot. Let us furthermore denote p(m) the 
probability density function (PDF) of the magnitude of 
the random sources. Then, the total average number of 
events (including the spontaneous sources and all their 
offsprings over all generations) during the time interval 
r with magnitudes larger than p is equal to 



where 



wt[<JZ)(aO+Q(aO] , 



(R)(fi) = / (R)(m, jj)p{m)dm 



(17) 



(18) 



is the average number of events of all generations with 
magnitudes larger than or equal to [i which are triggered 
by the spontaneous sources of all possible magnitudes 
above some lower threshold m s defined as the smallest 
magnitude of spontaneous events. In what follows, we 



assume that the spontaneous sources have their magni- 
tudes distributed according to a GR law 



p{m) = xe -x(™- m =) H(m - m s ) 



(19) 



with an exponent x possibly distinct from those of trig- 
gered events. The complementary cumulative distribu- 
tion function (CDF) of spontaneous sources magnitudes 
then reads 



(13) Q(jm) 



p(m)dm 



-x(p-m B ) 



[i > m s . (20) 



It is natural to introduce a magnitude threshold md 
of catalog completeness, i.e., only events with m > md 
are observed. Then, the total fraction of events above 
magnitude p, among all observable events in the time 
window r is given by the following normalization of (1171) 



F{n, m, d ) 



(R)(p) + Q(ti) 
{R)(m d ) + Q(m d ) 



(21) 



F(/i, m d ) can be interpreted as the complementary CDF 
of magnitudes of observable events. The corresponding 
PDF of observable events is then 



f(fi,m d ) = 



g(fi)+p(fi) 
(R)(m d ) + Q{m d ) 



9(fJ>) = - 



dfx 



(22) 

Expression l|17|) also allows us to obtain the fraction of 
triggered events 



n(/x) 



(#)(//) + QO0 



(23) 



which is nothing but the average branching ratio |27| . 

In the sequel, we apply relations (1221) and (|23fl to the 
standard ETAS model and to the self-similar Vere- 
Jones model ©■ 



IV. CRITIC ALITY CONDITION AND PHASE 
DIAGRAM 

A. General maximum eigenvalue condition 

Before going further, it is important to derive the con- 
ditions under which the branching process is not explo- 
sive, i.e., for which the stationary spontaneous sources 
point process generates a stationary sequence of trig- 
gered events. The condition derives in general from an 
eigenvalue problem (see |2^| and references therein). In 
the present case, it is known that equation 1|15H gives 
bounded stationary solutions if the largest eigenvalue p 
of the corresponding homogeneous equation 



plZ(m) 



X(m, m')lZ(m')dm' 



(24) 
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is smaller than 1 (p < 1). This corresponds to the sub- 
critical regime. The condition p = 1 defines the criti- 
cal regime and p > 1 gives the explosive super-critical 
regime. In our analysis of the eigenvalue problem of 
equation 1)24(1 , we shall restrict to physically meaningful 
eigenfunctions lZ(m) which are monotonically increas- 
ing functions growing no faster than the productivity 
law - e am : 



lim lZ(m)e 



< oo 



(25) 



An important point should be noted in relation with 
Vere- Jones's model. While in the standard ETAS 
model, the total number of triggered events above any 
arbitrary magnitude m& diverges at criticality, we will 
see that the average total number (R)(md) of observ- 
able events is finite in the critical regime of Vere- Jones' 
model, while the average total number (R)(rrid — > — oo) 
is itself infinite. Thus, while the process is critical when 
considering all events of any magnitude, it becomes sub- 
critical for events above a finite threshold m^. The hall- 
mark of such subcritical regime is that the fraction (|23[1 
for observable events is smaller than 1, i.e. n(p) < 1 for 
all p > irid- 



B. Criticality condition for the standard ETAS 
model 



For the standard ETAS model defined in section Hi Bl 
equation l|2"4"|) reduces to 



plZ(m) 



K(m')dm . (26) 



Let us look for a solution of this equation in the form 

K{m)=Ce Sm . (27) 

Substituting (|27() in 126fl shows that C ^ if and only 
if 5 = a while the corresponding eigenvalue is 



p = k / e( a -^ m dm 



n'(3 
13 -a 



(28) 



which is nothing but the average branching ratio n de- 
fined in (JSJl. This recovers the known fact that the 
condition p = n < 1 corresponds to the sub-critical 
regime associated with the solution 



(R)(m,p) = e «(™-mo)-/3(M-mo) (29) 

1 ~P 



of the non-homogeneous equation l)15[l. 

Substituting and (JTHJ) into JTHJ) and and 

assuming for simplicity that m s — mo, we obtain the 



distribution of the PDF of the magnitudes p of observ- 
able events 



(Rj e -@{rrtd-m a ) _|_ e -x(m d -rn a ) ' 



where 



(R) = (R)(m ) = 



1-p 



P = 

X - a 



(30) 



(31) 



is the average of the total number of aftershocks trig- 
gered by one spontaneous source of arbitrary magni- 
tude. In the particular case where the GR laws for 
the magnitudes of the spontaneous sources and of the 
first generation aftershocks are the same, i.e. if x = Pi 
then the PDF of the magnitudes of observable events 
given by (pTTTf> reduces to the pure GR law f(fi, m,j) = 
Pe-^ m - md) H{p - m d ). 
Substituting 



(R)(p) = 



,—PQt-mo) 



1-p 

into expression l(2"3")) yields 

P' 



GO*) 



n(p) 



p' + (1 - p)e(' 3 -x)(M-™ ) 



(32) 



(33) 



If, as before, x = P (p 1 = P)i then n{p) does not depend 
on p and is equal to the average branching ratio n = p. 



C. Criticality condition for the self-similar 
Vere-Jones model 

For the self-similar Vere-Jones model defined in sec- 
tion [^21 the homogeneous equation (|24|l reduces to 



p1Z(m) 



A(m - m')TZ(m')dm' , (34) 



where A(m) is given by expression Jfjjl. We search again 
a solution for TZ(m) of the form J27J, which yields the 
following eigenvalue p(S) as a function of 5 (shown in 
figure 1): 



p(S) = / \{m)e- 6m dm 



2nd 



d 2 -(5~ pf 



where 



1- (5-/3) 2 /d 2 ' 



2 k 



(35) 



(36) 



is going to play an important and physically intuitive 
role in the sequel. Note the major novelty compared 
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with the standard ETAS model: here, we obtain a con- 
tinuous spectrum of eigenvalues p(5) rather that the 
unique one l|28|) associated with S = a for the ETAS 
model. 

In expression (JSHJl, since 1 — (6 — (3) 2 jd 2 < 1, there is 
a solution with p{5) < 1 associated with the sub-critical 
and critical regimes only for s < 1. For s > 1, all 
eigenvalues p(S) are larger than 1, which corresponds to 
the explosive super-critical regime. We show below that 
the parameter s plays the role of an average branching 
ratio. 

For s < 1, the continuous spectrum of eigenvalues 
is indexed by 5 spanning the interval [f3 — h,f3\, for 
which p(S) < 1 characterizes the sub-critical and crit- 
ical regimes. The exponent h is determined by the con- 
dition p{[3 — h) = 1, whose geometrical determination 
is represented in figure 1. We rule out the possibility 
6 > [3, which leads to unphysical solutions as shown be- 
low. Given this spectrum of eigenvalues, the growth of 
lZ(m) given by (|27[) is controlled by the largest eigen- 
value p(/3 — h) = 1, when it exists, which leads to 



n{m) 



o {0-h)m 



dVT 



(37) 



We can now describe the phase diagram of the self- 
similar Vere- Jones model, shown in figure 2. 

• For s > 1, all eigenvalues p(8) are larger than 
1, corresponding to the explosive super-critical 
regime. 

• For s < 1 and d < j3, we accept (137(1 as a physi- 
cally appropriate eigenfunction only if it is mono- 
tonically increasing with respect to to, that is, if 
< h < (3. This leads to the condition 



d 



< K < 



d 



(38) 



In the range (|38[1 of parameters, there is always a 
solution, whatever the value < s < 1 of the form 
IpTTjl. associated with the unit eigenvalue. This 
is the critical regime, which is associated with a 
finite range of parameters < s < 1 and < h < 
f3 corresponding to the domain l(38|l . indicated in 
figure 2. 



• For 



Q < K < - ( d 



d 



(39) 



which is only possible if d > /3, we obtain h > (3, 
that is, S < 0, which corresponds to a seismic ac- 
tivity which is a decreasing function of the mother 
magnitude. This is the hallmark of the subcrit- 
ical regime. We will see below that the average 
of the total number of offsprings is finite in this 
sub-critical regime. 



Figure 2 summarizes this phase diagram and delineates 
the domains of existence of the sub-critical, critical and 
super-critical regimes in the plane of parameters (k, d) 
for a given [3. 

The physical meaning of this classification is ob- 
tained by examining the equation for the average num- 
ber (R) (to, p) of aftershocks which are triggered by some 
mainshock with magnitude m. This equation is expres- 
sion H15fl written for the self-similar Vere- Jones model: 



(R)(m,p)= / \(m-m')(R)(m',p)dm / +(R 1 )(m,p), 

J — OO 

(40) 

where 

P CO 

(R 1 )(m,p)= / \(m-m')dm! (41) 

is the average of the number of corresponding first- 
generation aftershocks. Let us introduce the auxiliary 
function 



S(m, p) = — 



d(R)(m,p) 
dp 



(42) 



Using (jl^l) . S(m,p) is solution of 

/>OC 

S(m,p)= / X(m—m')S(m',p)dm'+X(m—p). (43) 



For its structure, it is clear that the solution of this 
equation depends only on the difference between to and 
p, so that we write 



S^to, p) = S(m — p) 



(44) 



The value of (R) (to, p) is obtained from S*(m, p) by us- 
ing 



(R)(m,p)= / S{m,p')dp' = / S(x)dx . (45) 
Let us solve the equation 

S(m) = / A(m - m')S(m')dm' + A(m) (46) 



for the auxiliary function S(m). Applying the two-sided 
Laplace transform to this equation yields an equation for 
the Laplace transform 



S(u) 



S(m)e- um dm 



which takes the following form: 

A(u) 



S(u) 



1 - X(u) 



(47) 



(48) 
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The Laplace transform X(u) of the kernel A (to) given by 
© can be explicitly calculated as 



X(u) 



2nd 



d 2 -(u- [3) 2 ■ 
Substituting expression lUD into gSJ yields 
2nd 



S(u) 



h 2 - (u-f3) 2 ' 
Its inverse Laplace transform 



h l =d 2 - 2nd 



s <"*> - h 



&-\-ioo 



S(u)e um du 



(3—ioo 



gives finally 



S( m ) = — e 0m-h\ m \ 
h 



(49) 



(50) 



(51) 



(52) 



We can now use (|52() in (|45() to characterize the aver- 
age of the number of aftershocks in the sub-critical i|39|) , 
critical lj38J) and super-critical (s > 1) regimes. 



• In the subcritical case h > /3, we can take fi 
and still obtain a finite limit 



(R) 



lim (R)(m, //) 

[A — > — OO 



S(x)dx 



2nd 

h 2 -(3< 



< OO 

(53) 

for the average of the total number of aftershocks 
of all generations triggered by an arbitrary spon- 
taneous source. Thus, in the subcritical regime, 
the averages of total number of both observable 
and unobservable events are finite. 

In the critical regime (|38[l , we find that the average 
(R) (to, /i) of the total number of events above any 
finite magnitude threshold /i is finite. In contrast, 
the average of the total number of aftershocks 
(including the unobservable tiny events of mag- 
nitudes m! — > — oo) becomes infinite: (R) = oo. 
It is remarkable that we have at the same time 
(R)(m,iJ,) finite for any /i > — oo and (R) = oo. 
In real data, we only observe (R)(m, /i). The un- 
derlying criticality is thus unobservable: due to 
the special cancellations in the self-similar Vere- 
Jones model, the infinite swarm of tiny events 
form an unobservable sea of activity, whose ob- 
servable consequences lie in the finite activity at 
finite magnitudes. Thus, in a sense, this critical 
regime can actually be decomposed into an 
effective sub-critical regime for s < 1 and a criti- 
cal point reached at s = 1 for observable events. 
Indeed, as s — > 1 (i.e. if h — * 0), expressions (|45fl 
and il~>2[| for the average number of observable af- 
tershocks triggered by an arbitrary spontaneous 



source tends to infinity. This shows that the crit- 
ical regime for observable events corresponds to 
s = 1 and confirms the interpretation of s as the 
effective branching ratio for observable events. 

As s — > 1 (h — > 0), (R)(m,fi) tends to infinity for 
any to and fi, confirming that the super-critical 
regime corresponds to s > 1. 



D. Criticality condition for the generalized 
Vere-Jones ETAS model 

We now analyze the critical conditions for the gen- 
eralized Vere-Jones ETAS model defined by expression 
(0 in section III Al It is convenient for the analysis to 
represent the intensity A(to,to') in J5J) by 



A(to,to') = e Q(m ~ m 'MTO,TO') 

where 

v(m,m!) = Ke ^-PW-d\ m - m '\ 
Let us introduce the auxiliary function 
U(m) = K(m)e~ am , 



(54) 
(55) 
(56) 



and rewrite the homogeneous equation l|24|l in the form 



pU(m) 



U{m)v(rn ) m!)drri . 



(57) 



From the definition of v{m, to') in l|55|l . we see that the 
following equality holds: 

d 2 is(m m>) = _ e{a _ 0)m _ ^ 

dm* 

(58) 

where 5(x) is the Dirac-delta function. By differentiat- 
ing equation l|57l) twice with respect to to, the identity 
l|58|) leads to the condition that, if U(m) is a solution 
of Eq. (|57(l . then it should also satisfy the differential 
equation 



d 2 U{m) 



dm 2 



pd 2 - 2Kde {a ~ 0)rn 



U(m) , 



to > too . 
(59) 

Thus, to determine the eigenvalue of equation (|57|l . our 
strategy is to search for a solution of equation l|59|l . 
which at the same time satisfies the integral condition 

/•OO 

pU(m ) = K / U{m)e [a - p)m ~ d{m - ma) dm . (60) 



which is derived from l|57|l . It is straightforward to check 
that this strategy leads to solving l|57|l . In addition, we 
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need to impose that the eigenfunction lZ(m) is monoton- 
ically increasing as a function of m such that the condi- 
tion (|25[1 is satisfied. In terms of the auxiliary function 
U(m), this implies that 



lim hOrn) < oo 



(61) 



Thus, our problem is to find a solution of (|59|) with the 
two conditions (ffitHl and ljfil|) . 

Before using this formalism for the generalized Vere- 
Jones ETAS model, it is instructive to see how it per- 
forms on the standard ETAS model (corresponding to 
d = 0). For d = 0, equation (|S*9")l reduces to 



d 2 U{m) _ 
^ dm 2 



(62) 



Its solution satisfying condition (|61|l is U = C — const. 
Substituting it into the integral relation 1)61) |) recovers 
the known expression for the unique eigenvalue equal to 
the average branching ratio: 



P : 



e (a-/3)m Q 

/3-a 



1 i 

K 



7-1 



/3 

7 = - 
a 



(63) 



The critical regime corresponds to the set of parameters 
obeying 



K' = fC(a,j3) = 1 



(64) 



while the domain of subcritical regime corresponds to 
k < JC(a, (3). Similarly, to explore the condition for crit- 
ically in the generalized Vere-Jones ETAS model, we 
just need to put p = 1 in (|59|l and search for the func- 
tion 



K = JC(a, (3, d) 
such that the homogeneous equation 



(65) 



d 2 U(m) 
dm 2 



d 2 - 2Kde {a ~ p)m 



U{m) 



m > mo 



( 66 ) 

has a nontrivial solution increasing with m, which sat- 
isfies the integral equality (|60|l expressed for p = 1 . The 
sub-critical regime then corresponds to the set of pa- 
rameters such that k < K,(a,0, d). 

In order to solve H66fl for the critical case p = 1, we 
introduce the new function v(y) such that 



with 



V = Qe 



-(m—mo)d/e 



U(m) = v(y) , 

VMdp 

g = ~n 

p — a 



(67) 



2d 



(68) 



Note that the parameter e quantifies the transition from 
the ETAS model (obtained for e = 0) to the self-similar 
Vere-Jones model (obtained for e — > +oo). For small 
(resp. large) e, the generalized Vere-Jones ETAS model 
is close to the standard ETAS (resp. self-similar Vere- 
Jones) model. Equation (|66|1 for U (m) translates into 
the following Bessel equation for v(y): 



id 2 v{y) dv(y) 



dy 2 



dy 



+ (y 2 - e 2 )v(y) = 



(69) 



It follows from l|61|l and from the definition (|68|l of y 
that the solution of (|69(l has to satisfy 



lim v(y) < oo 



(70) 



The integral condition l|60(l imposes in addition that 

v{Q) = \e- 1 J\(y)y* +1 dy . (71) 
The solution of Eq. l|69fl has the form of a Bessel function 



v{y) = AJ e {y) 



(72) 



where A is a constant. Substituting 172(1 into (|71|l and 
using the known recursion relation 



y< +1 J e (y)dy = y^ 1 J e+1 (y) 



(73) 



between Bessel functions, we obtain the implicit equa- 
tion for the variable g 



(74) 



which determines the set of parameters corresponding 
to criticality (p — 1). 

Further insight can be obtained by determining the 
leading behavior of g for small e (quasi-ETAS model) 
and large e (quasi-self-similar Vere-Jones model). For 
this, we introduce the new auxiliary parameter 



J, 



2k'/3 



and rewrite equation l|74|) in the form 

J e + l(«/>) 

For e<Cl, we use the expansion 



(75) 



(76) 



J e+ i(x) x 



2 + e 



(x « 1) (77) 
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which gives the solution 

V * j- (2 + 6) (6 « 1) 



(78) 



In term of the original parameters k, a and d, this gives 



a d 



(79) 



as the relation expression the critical regime p = 1 of the 
generalized Vere- Jones ETAS model. Expression (|79|l 
differs from its counterpart (|64() obtained for the stan- 
dard ETAS model by the correction term d/P, which 
describes a kind of broadening of the subcritical regime 
due to the magnitude localization effect (d > 0). 

In the other limit e » 1 corresponding to the quasi- 
Vere- Jones model, the solution of the criticality condi- 
tion (|76[l can also be obtained asymptotically by expand- 
ing the Bessel functions in the neighborhood of their first 
zero, i.e., by searching for ep close to v = v(e) defined 
by J e [f(e)] = 0. The corresponding expansions in Taylor 
series of the Bessel functions are 

J e (x) ~ J' e (v)(x-v) , J e +i(x) ~ J e+ i(v)+J' e+1 (v){x-v) 

(80) 

It is well-known from the theory of Bessel functions that 



J' e (v) 



-Je+lM 



1 



(81) 

Thus, for e > 1, expressions (I8UII and H81|) allow us to 
replace the criticality condition (|76|l by its asymptotic 
expression 

2v(v — eip) 



v+ (e + l)(v- eip) 
Solving this equation in ip yields 
v(e) 



4> 



(82) 



(83) 



Using the known asymptotic relation 

u(e) ~ e + 1.856 e 1/3 + 1.033e- 1/3 (84) 

for the first-zero of the Bessel function of order e, the 
criticality condition formulated in terms of the original 
parameters k, a, d then reads 

-'-^ + ^ 1/3 (/3-«) 2/3 - (85) 

The first leading term n' « ^ recovers the critical con- 
dition s = 2k/ d = 1 for the self-similar Vere- Jones 
model. The last term in the r.h.s. of l|85l) thus provides 
the first correction to the exactly self-similar Vere- Jones 
model when a ^ (3. Figure 3 plots the numerical solu- 
tion of equation l|76[l as a function of e together with its 
two asymptotics regimes (|78fl and (|83|) with l|84|) . 



V. DISTRIBUTION OF EARTHQUAKE 
MAGNITUDES FOR VERE-JONES'S 
SELF-SIMILAR MODEL 

The Gutenberg-Richter distribution of earthquake 
magnitudes is perhaps the most ubiquitous and docu- 
mented statistical property of earthquake catalogs. It 
is thus of great interest to investigate the prediction of 
Vere- Jones self-similar model of this quantity. We an- 
alyze in turn two versions of this statistics, first when 
the magnitude of the mainshock is known and second 
when summing over all possible magnitudes as in a large 
catalog. Both forms of the distributions have been in- 
vestigated in the empirical literature. 



A. Distribution of the magnitude of aftershocks 
triggered by a mainshock of given magnitude rn 

Reasoning as in section fill Bl shows that the PDF's 
of the magnitudes p of first-generation aftershocks and 
of all aftershocks of all generations of a mainshock of 
magnitude m are proportional respectively to Si(m, p) 
and S(m,p) given by 

d{R x ){m,p) d(R)(m,p) 
Di(m, p) = , b(m, p) = - 



Op 

Their expressions are 



dp, 



(86) 



(87) 



Si(m,p) = Ke -^-m)-d\n-m\ ^ 
S (m,p) = ^ e -^-m)-h\„-m\ 

h 



Thus, the two branches p < m with exponent (3 — d 
and p > m with exponent (3 + d of the PDF S\{m,p) 
of first-generation aftershocks are renormalized in two 
other branches of the same form, with d renormalized 
into h = d\/l — s. The main point is that only S(m, p) 
is observed in real catalogs since one does not have the 
luxury of being able to distinguish between the different 
generations of triggered aftershocks. Expression (|88|) 
shows that the PDF S(m, p) may be quite close to a 
single power law, even if the difference between the two 
branches of Si(m, p) is significant, as long as the critical 
parameter s is close to 1 (for s — > 1, S(m, p) ~ e _ ^ M and 
is a pure power law). In other words, an observable PDF 
of aftershock magnitudes close to a single power law is 
compatible with strong deviations from a pure power law 
for first-generation events, due to the renormalization 
effect over all the generations which effectively mixes 
up the two branches sufficiently close to criticality s — + 
1. Thus, as in previous studies of the standard ETAS 
model 11(1 MJ. comparison between the data and these 
predictions suggests that the earth is operating close to 
a critical point. 
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B. Distribution of earthquake magnitudes over all 
events 

The Gutenberg-Richter law for the distribution of 
the magnitudes of earthquakes is generally a statistical 
property established for a large space-time-magnitudc 
domain, without restrictions. It is thus interesting to 
ask what Vere- Jones' self-similar model predicts for the 
distribution of magnitude of a large stationary sequence 
of events triggered by a steady-state influx of sponta- 
neous sources. The answer is given by expression (1221 
for f(/j,rrid)- We thus need to make explicit the depen- 
dence of /(/x, ma) on the parameters of the model. This 
problem depends on four key parameters: x is the ex- 
ponent for spontaneous sources defined in (|19[) : [3 is a 
crucial parameter in the definition of the PDF of mag- 
nitudes of first-generation aftershocks; h and s appear 
in the condition for criticality. 

Let us first investigate the contribution g([i) of the 
events triggered by the spontaneous source. For Vere- 
Jones' self-similar model, we have 



S(m — [i)p(m)dm 



(89) 



Substituting (JTHJ) and in 1^ yields 

g{ l x)=p{lJL)K(m s - l x,p-x) , (90) 

where K(m s — /j.,/3 — \) describes the deviation of the 
PDF of the magnitudes of triggered events from the 
PDF of the magnitudes of the spontaneous sources. This 
function is given by 



K{y,z) = ^ I" e zx - h ^dx 
h J y 



Its explicit expression is 
2nd 



K(y,z) = 



Kd 



Jz+h)y 



h? - z 2 h(h + z) 



(91) 



(92) 



The corresponding complementary cumulative distribu- 
tion, i.e., the total number of triggered events with mag- 
nitude larger than /i, is equal to 

/>oo />oo 

(R)(n) = / g(m)dm = / p(m)K(m s —m, (3—x)dm . 
J n J i± 

(93) 

We need to distinguish three cases as g{^) given by 
(I9L)|) is qualitatively different for different values x> P an d 
h. 

• For x < — h, g(fJ-) = oo for any /i. This re- 
sults from the fact that the spontaneous sources 
with large magnitudes dominate the production 
of triggered events in this case. This super-critical 



regime can be tamed with the introduction of an 
upper magnitude cut-off m mayi but is not investi- 
gated further here. The limiting case x = P ~ h 
gives a criticality condition for the observable 
events in the framework of the self-similar Vere- 
Jones model. 



Consider the regime 



l3-h<x<f3 + h, 



(94) 



and the limit m s — > — oo corresponding to y — > 
— oo, for which K(y, z) has the following asymp- 
totic dependence 



K(y,z) 



2nd 



h 2 



y 



-oo . 



(95) 



The corresponding asymptotic expressions for the 
PDF g(ii) given by l|9Q(l and the complementary 
cumulative distribution given by l|93|l are 



2nd 



h 2 



2nd 



h 2 



■Q(fi) . (96) 



Substituting lO and ^ into yields 

the distribution of the magnitudes of all observable 
events (in the limit m s — > — oo) 



f(fi,m d ) =xe~ x{ " d) (V>m d ) 



(97) 



Remarkably, in this regime (|94J) . the observed 
Gutenberg-Richter distribution is predicted to re- 
veal only the exponent of the distribution of the 
spontaneous sources and to be blind to the expo- 
nents P, d of the distribution of triggered events. 
In other words, the PDF of the magnitudes of 
all observed events reproduces that of the spon- 
taneous sources given by ljT§|) . 

Substituting (|96|l into (|23|l yields the fraction of 
triggered events whose magnitudes are large than 
fj,, among all analogous events, 



n(n) = n = 



2nd 



d 2 -{P~ x) 2 



(98) 



which is found independent of fi. If the same 
distribution of magnitudes describes the sponta- 
neous sources and the triggered events (/3 = x)> 
then n = s. Thus, the critical regime corresponds 
to s = 1, confirming the interpretation of s as 
equivalent to the critical branching ratio of Vere- 
Jones' self-similar model. However, if /3 7^ x an d 
|/3 — x| — > h, then we obtain n — > 1 even for s < 1 
(h > 0): essentially all the events are triggered. 



• Consider the regime 



X > P + h 



(99) 
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In this case, the function K is closely approxi- 
mated by its asymptotic behavior 



K(y,z) 



K(l 



h(h + z) 



(100) 



The corresponding asymptotic behaviors of g{n) 
given by (|§U|l and {R){p) given by for large 
(i — m s are 



and 



-(/3+7i)(^-m s 



~ h( X -(3- h) 

i^tx 

h( x -f3-h)((3 + h)' 



-(f3+h)(p- mB ) 



(101) 



(102) 



Substituting these two expression into 122fl and 
(EH) yields 

f(n, m d ) ~ (/? + h) e -(0+h)(n-™ d) , „ ~ i j (103) 
asymptotically for /i — m s — > +oo. 

VI. SUMMARY AND DISCUSSION 

Our results can be summarized as follow. 

• We have clarified and quantified the conditions un- 
der which the self-similar Vere- Jones model as well as 
a more general version (which contains both the stan- 
dard ETAS and Vere-Jones version as special cases) are 
critical, subcritical and supercritical. Only the subcrit- 
ical and critical regimes give a stationary process in the 
presence of a non-zero flux of immigrants. 

• We have shown that the concept of an average 
branching ratio s, defined as the average number of 
daughters of first generation per mother of magnitudes 
above a finite magnitude threshold) holds for Vere-Jones 
model in a broad domain of parameters. Remarkably, s 
is found independent of the magnitude threshold used. 
However, the average branching ratio loses its mean- 
ing when the magnitude threshold is pushed to — oo (in 
other words, when it is removed), as the existence of 
arbitrary small events allowed in this model dominates 
and makes the average divergent. Since empirical cata- 
logs are always characterized by a minimum magnitude 
md of completeness, our results apply directly and show 
that it is possible to have an infinite number of - unob- 
servable but still important for the cascade of trigger- 
ing - events per mother together with a finite average 
branching ratio for observable events. 

• Vere-Jones' model is defined by the Gutenberg- 
Richter (GR) magnitude distribution for first-generation 
events triggered by a source of magnitude m having two 
branches: for aftershocks magnitudes ml < m, the GR 



exponent is (3 — d, while it is /3 + d for ml > to. We 
have shown that, accounting for the contributions of all 
generations of triggered events, this GR distribution is 
renormalized into another two-branches law: for after- 
shocks magnitudes ml < to, the renormalized GR ex- 
ponent is (3 — h, while it is (3 + h for m! > m, where 
0<h= d^/T~s < d for < s < 1. 

• Since only the renormalized GR is observable in the 
analysis of empirical catalogs, the fact that the literature 
reports a single exponent (3 for the GR distribution im- 
plies that, assuming that Vere-Jones model is a correct 
description, h is small so that the difference between the 
exponents (3 — h and (3 + h is within the empirical un- 
certainties and variations from catalog to catalog. This 
implies that either d is small or s is close to 1 (condi- 
tion describing the boundary between the subcritical to 
supercritical regimes of the Vere-Jones model) or both. 
We note that the prediction that two distinct exponents 
h and d characterize respectively the GR law over all 
generations and the GR law over the first generation 
of triggered events could be in principle tested empiri- 
cally using an adaptation of the statistical declustering 
method developed by Zhuang et al. 0, El to Vere- 
Jones's model. In this respect, Zhuang et al. have 
already found that the magnitude distribution of the 
triggered event depends on the magnitude of its direct 
ancestor, with an exponent smaller for large events (in 
contradiction with other less sophisticated studies us- 
ing more arbitrary space-time windows |l8t Il9l ] ) : this is 
roughly consistent with our prediction with Vere-Jones 
model, as the observable distribution for large (resp. 
small) ancestors is weighted more by the m! < m (resp. 
m! > to) regime associated with exponent (3 — h (resp. 
13 + h). 

• The distribution of magnitudes over a stationary 
catalog (obtained by summing over an average steady 
flow of spontaneous sources) is found universal (inde- 
pendent of 0, d and the other parameters) and a pure 
GR with exponent equal to the exponent x °f t ne spon- 
taneous sources (which in full generality is allowed to 
be different from the exponent (3 involved in the dis- 
tribution of triggered events) in a large domain of the 
parameter space. This implies that, if the exponent of 
the spontaneous sources is different from the exponent 
of triggered events, the physics of cascades of triggering 
in the self-similar Vere-Jones model implies that only 
the former exponent is observable in global catalogs. 
Ag ain, the statistical declustering method of Zhuang et 
al. [16l ri7j | should be able to test this prediction. In this 
respect, we note that Zhuang et al. find that the back- 
ground events have a larger exponent than the triggered 
events: x > P- 
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(3-d p-h p p+d 

Fig. 1: Dependence of the eigenvalue p given 
by l|35fl associated with the eigenfunction lZ(m) 
1271 1 . as a function of S for the self-similar Vere- 
Jones model. 




p 

Fig. 2: Phase diagram in the plane of param- 
eters (k, d) at fixed (3, for the self-similar Vere- 
Jones model, showing the domains of existence 
of the sub-critical l|39fl . critical l|M8fl and super- 
critical (s > 1) regimes. 
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Fig. 3: Numerical solution of equation 1761 as 
a function of e defined in l|68|l compared with 
its two asymptotics regimes 178H and 1831 with 
iJHU, giving the condition for criticality p — 1 
for the generalized Vere- Jones ETAS model. 



8 



